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TECHNICAL NOTE 3343 


SUBSONIC EDGES IN THIN-WING AND SLENDER-BODY THEORY 


By Milton D. Van Dyke 
SUMMARY 


A simple technique is presented for correcting the results of thin- 
wing theory near round or sharp edges at which the normal component of 
free~stream velocity is subsonic. The flow in planes normal to such 
edges is actually brought to rest, but thin-wing theory instead predicts 
infinite velocities. Furthermore, corresponding to the circulatory flow 
around leading edges, thin-wing theory predicts additional singularities 
that do not actually exist for round edges and exceed their true strength 
for sharp edges. All these singularities grow worse if one formally 
attempts to improve the solution by iteration. To correct these defects, 
the formal thin-wing solution is here rendered uniformly valid by compar- 
ing with exact solutions for simple profiles having the same nose shape. 

In this way Lighthill's rule for correcting the surface speed on round- 
nosed airfoils in incompressible flow is extended to higher approximations, 
to compressible flow, to three-dimensional wings, and to sharp edges. 
Corresponding results for slender bodies of revolution are considered 
briefly; in particular, the rules for round-nosed airfoils in incompressi- | 
ble flow are shown to apply also to round-nosed bodies of revolution. = 


INTRODUCTION 


The Linearized theory of thin wings, which has proved so fruitful 
for both subsonic and supersonic flows, is known to fail near leading and 
trailing edges if the component of free-stream velocity normal to the edge 
is subsonic. The flow in a plane normal to such an edge is actually 
brought to rest at some point near the edge, but linearized theory predicts 
infinite velocities instead.t Furthermore, there is flow around the lead- 
ing edge except at one angle of attack, and linearized theory then predicts 
an additional singularity that does not actually exist for round edges 
and exceeds its true strength for sharp edges. 


These local failures are sometimes of little practical consequence, 
Since the singularities are integrable and the lift is determined correct 
to a first approximation. However, it has long been known that special 
a a A a a r E AA 


*At round edges, linearized theory predicts only finite speeds on the 
surface, but the velocity of the stream approaching the stagnation point 
becomes infinite. 
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attention must be paid to the singularity in calculating the ärag of thin 
lifting wings (ref. 1). More recently, R. T. Jones has pointed out that 
the same is true for nonlifting wings and has shown how to account for 
the singularity in evaluating drag (ref. 2). In any case, the calculated 
velocities and pressures are incorrect near such an edge. 


If one attempts to refine the theory by formally iterating upon the 
first approximation, the local failure leads to serious difficulties. 
In incompressible flow, the singularities are found to increase in 
strength at each stage of the iteration, so that the solution, though 
improved elsewhere, becomes increasingly inaccurate near the edge. More- 
over, for round edges these higher terms give nonintegrable contributions 
to the aerodynamic forces. Even more serious difficulties may arise in 
compressible flow, where the infection spreads in some cases, so that the 
formal second approximation is incorrect not only at the edges but over 
the entire airfoil surface. 


The remedy for all these difficulties is clearly to base the itera- 
tion upon a first approximation that is valid.even near the edges. This 
would appear to require abandoning at the outset the assumption of small 
disturbances, which so greatly simplifies analysis. Fortunately, however, 
Lighthill has found (ref. 3) that for round edges in incompressible flow 
all the results of small-disturbance theory can be salvaged. He shows 
that a very simple rule converts the formal second-order solution of thin-~ 
airfoil theory into an approximation that is uniformly valid near the 
edge. 


Lighthill's result follows from application of his general technique 
for rendering approximate solutions to physical problems uniformly valid 
(ref. 4). This consists in determining progressively with each approxi- 
mation a straining of the coordinates that ensures uniform convergence. 
For a round-nosed airfoil in incompressible flow the straining is, to the 
first approximation, a contraction of abscissas that shifts the leading 
edge by half its radius, as was noted by Munk in 1922 (ref. 5). The uni~ 
formly valid solution is then the same function of the strained coordi- 
nates as the thin-airfoil solution was of the unstrained coordinates. 


It is not yet clear how far Lighthill's technique can be extended. 
Fox has indicated (ref. 6) that it fails for higher approximations, other 
nose shapes, or compressible flows, but in a recent series of seminars 
at California Institute of Technology, Tsien suggested that proper appli- 
cation of the technique will lead to success.: In any case, a quite dif- 
ferent technique is used here, which is particularly suited to the study 
of edges, and for that special problem is simpler than Lighthill's general 
procedure. It consists in comparing the exact and thin-wing solutions 
for simple shapes that approximate the airfoil_in the. vicinity of the 
edge. Thus Lighthill's rule for round noses is reproduced by considering 
incompressible flow past a semi-infinite parabola. This comparison tech- 
nique is then readily extended to higher approximations, compressible | 
flows, three-dimensional wings, sharp edges, and slender bodies of revo- 
lution. Attention is confined here to the practical matter of calculating 
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surface speeds, but the whole flow field can be corrected by the same 
method. The technique is applied here only to edges that are round 
(finite radius of curvature) or sharp (finite slope), since these occur 

on most existing profiles. However, the procedure can, in principle, be 
extended to other edge shapes and to slender bodies of other than circular 
cross section. 


The analysis of this paper is intuitive in spirit. Where convenient, 
some of the results have been deduced by appeal to physical reasoning. 
It is presumed that at these points the analysis could be bolstered by 
more rigorous mathematical arguments. 


The author is indebted to R. T. Jones for many helpful discussions 
throughout the course of this work. 


ROUND-NOSED ATRFOILS IN INCOMPRESSIBLE FLOW 


The technique to be used can be introduced by restricting attention 
first to an uncambered two-dimensional airfoil at zero angle of attack in 
incompressible flow. Thin-airfoil theory breaks down in the vicinity of 
the leading edge. There, a round-nosed airfoil can be approximated by a 
parabola. Accordingly, consider the exact velocity distribution over a 
parabola of leading-edge radius P in a uniform stream parallel to its 
axis (sketch (a)). The flow speed on 
the surface can be found by conformal Y 
mapping, or simply by Munk's rule t 
(ref. 7) that the surface velocity on i 
any ellipsoid subjected to uniform 
flow along a principal axis is the Y 
projection of the maximum velocity 
onto the plane tangent to the surface. 

For a parabola the maximum is the free- 
stream velocity, so that with x 
measured downstream from the nose the Sketch (a) 


surface speed ratio is 
1/2 
= (- a, ) (1) 
+ p/2 


(All symbols are defined in the Appendix.) Expanding this expression for 
small p/x yields the formal series 


YutV2Z ox 






cya 


LE S F i 
q p 
U = 1 — lx + e o o (2) 
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and this must be the series that would be given by thin-airfoil theory. 
(It is actually the second-order thin-airfoil solution - the first-order 
term being zero for a parabola - since the nose radius of any profile 
varies with the square of its thickness.) It is clear how the spurious 
Singularities of thin-airfoil theory arise at the leading edge. 


The formal second-order thin-airfoil solution "q5 for a parabola 


can be converted into an approximation q that is uniformly valid by 
multiplying it by the ratio of the exact result to the series expansion: 


1/2 
(5) | gn 
x + p/2 2 


U 


(3a) 


GOs 


ec O 


hx 


Simplifying this insofar as possible, retaining only terms of order o/x, 
gives the rule 


q x ae 1 “dp p 
eae pL 3b) 
= oe o U ix l 


It can now be argued that this rule, derived for a parabola, applies 
with some degree of approximation to any uncambered airfoil of finite 
leading-edge radius Pf at zero angle of attack. Near the leading edge, 
where the flow disturbances are large, the exact speed on any airfoil is 
nearly that on a parabola of the same nose radius; hence, the rule pro- 
vides a first approximation to the disturbances there. Far from the edge, 
the disturbances are small and, in general, different from those for a 
parabola; there the ratio in equation (3a) is essentially unity, so that 
the rule reproduces_the second-order solution for the airfoil, as it 
should. Hence, the rule yields an approximation that is uniformly valid. 
(It is not, of course, valid near a second edge; combined edges will be 
considered later.) 


The rule clearly applies also to 
cambered airfoils at the “ideal” angle 
Ler of attack for which the stagnation 
point coincides with the vertex. How- 
ever, at other angles of attack the 


* flow field includes a component of 
circulatory flow around the edge 
— (sketch (b)). Because of the linear 
q nature of incompressible flow, this 


a, component can be treated separately. 


The amount of this component is pro- 
Sketch (b) portional to the sine of the angle of 
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attack ao measured from the ideal 
angle. The thin-airfoil series for 
this component can be corrected by a 
comparison with the circulatory flow 

past a parabola (sketch (c)). Con- 

formal mapping or Munk's rule? gives 

the speed on the upper and lower 

surfaces of the parabola (counted 

positive away from the vertex) as 


1/2 





GI 


= tk sin ao ( l 


EET (4) Sketch (c) 


Here the factor k is a coefficient of order unity that relates the cir- 
eulatory flows on the parabola and the actual airfoil; it depends on the 
entire airfoil shape (and the trailing-edge condition), but its value is 
not required here. Expanding this expression for small o/x yields the 
formal thin-airfoil series for the circulatory component, which is, to 
second order in thickness and angle of attack, 


n" q> kao 


Me E (5) 
U y xX 
As before, the ratio of the exact result to its series expansion is a 


factor that serves to correct the thin-airfoil solution for any round- 
nosed airfoil. Hence the rule for the circulatory component is 


oa oF 9 


Now it happens that this is essentially the same as the rule for the 
other component (eq. (3b)). Consequently, the two rules can be combined. 
It follows that the rule of equation (3b) applies to a round-nosed air- 
foil at any angle of attack. It is, in fact, precisely Lighthill's rule, 
which he found using his technique for rendering approximate solutions 
uniformly valid. 





Gha 
l 


The accuracy of this rule near the nose of an airfoil depends upon 
how closely the nose is approximated by a parabola. Lighthill has 


“Munk's rule can be used here by recognizing that the circulatory 
flow about a parabola is identical with the limit of the crossflow past 
an ellipse as its thickness ratio vanishes. 
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considered airfoils described by analytic arcs; which means that near the 
leading edge ah 


y= e(+F,x7/* + Cx $F, x°/* teso‘) (7) 


€ being a small parameter and the Fp and Cy coefficients of order unity 
associated with thickness and camber, respectively; also, the angle of 
attack is assumed to be of order €. Now the example of the parabola 
shows that the problems of nonuniformity are confined to, a region propor- 
tional to the leading-edge radius, that is, in which x/ is O(e#). In 
this region the analytic airfoils of equation (7) coincide with parabolas 
except for terms of O(x?/#), which is O(¢€). It seems reasonable to assert 
that the surface speed will be approximated to the same order of accuracy 
as is the airfoil shape. This means that applying the rule to the second- 
order thin-airfoil solution yields a uniformly valid solution that is in 
error by O(€) near the leading edge (though of course only by O(e€8) far 
away). This agrees with Lighthill's conclusion, based on more rigorous 
analysis. However, initially uncambered airfoils (for which Ci = 0) are 
parabolic except for terms of O(x) = O(e#), so that in this case the 
solution rendered uniformly valid is in error by at most O(e*) near the 
nose. 


It would at first sight seem appropriate to say that Lighthill's 
rule yields a solution for analytic airfoils correct only to zero order 
near the nose (first order in the absence of camber), since only terms 
of O(1) are found correctly, and this is the order of the undisturbed 
stream. However, this occurs where the perturbations themselves are 
O(1) - for example, the surface speed one leading-edge radius from the _ 
nose is 0.707 of the free stream ~ so that the leading term in the veloc- 
ity disturbance (and hence inthe pressure coefficient) is given correctly. 
With this in mind, we shall say that Lighthill's rule renders the thin- 
airfoil solution valid to first order near the nose for cambered airfoils 


and to second order for initially uncambered ones. 


If only first-order terms are retained, Lighthill's rule becomes 


simply 
U X + p/ U. 


This converts the linearized solution into a uniformly valid first 
approximation. Results differing from this only by higher-order terms 
were given earlier by Goldstein (Approximation II of ref. 8) and by 
Riegels (ref. 9). It will be seen later that Riegels! form happens to be 
both simpler and more accurate than Lighthill's. 
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It seems desirable to refine the rule so that it is correct to 
second order ~ that is, to O(€2) - even near the nose of a cambered 
airfoil. This is readily done by 
considering flow past a "“cambered” 
or sheared parabola (sketch (d)) given 
by the first two terms of equation (7). J x 
Shearing transforms a parabola into \ 
another parabola whose vertex and nose f 
radius are changed only to higher \ x 
order. Hence the previous rule of RÉ i 
equation (3b) applies if x is isc [ 
replaced by the oblique abscissa X. 
Now it is easily shown that for a 
point on the surface of the parabola Sketch (d) 


X = x +eC,,/20x +... (9) 


where €Cı is the initial slope of the camber line. Hence the rule 
becomes 





do x t eCi/2px VF / "a 
== (10) 
x + 0/2 + €C,,/20x- 


Where, again, the upper and lower’ signs pertain to the upper and lower 
Surfaces of the airfoil. ; 


This rule, which renders the thin-airfoil solution uniformly valid 
to second order, is probably all that will ever be required in practice. 
However, it can easily be extended still further if necessary. For 
example, consider an uncambered analytic airfoil. Its nose can be 
approximated by an ellipse except for terms of 0(x°5/@), which is O(e>) 
in the region of interest. Hence, the following rule, derived from con- 
sideration of the exact flow past an ellipse, renders the thin-airfoil 
solution uniformly valid to fourth order: 


ra A 


ig 


LESIE i; 2 
A Sg l2 -É (Eo ) 
-== (=> T - BA ) 9 A= T ( x + 6FoFı (11) 


The fourth-order result for cambered airfoils can be obtained from this 
using the fact that a sheared ellipse is an ellipse; this will be given 
later when alternative forms are considered. 


SAn ellipse corresponds to F,<0; application of this rule to cases 
for which F,>0 rests on an appeal to the principle of analytic continu- 
ation. 
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From these rules, corresponding rules follow at once for any of the 
velocity components. For example, Lighthill's rule for the horizontal 
velocity component at the surface can be deduced directly from equa- 
tion 3(b): 


: “cme (H+ = £) te) 


In calculating pressures from velocities, the full Bernoulli equation 
must be used since the disturbances are not small. 


REMARKS ON THE METHOD 


The essence of the technique used here is to consider the exact 
solution for some simple shape that approximates the airfoil in the 
vicinity of its edge, where thin-airfoil theory breaks down. The formal 
thin-airfoil solution can then be corrected by multiplying it by the 
ratio of this locally exact solution to its series expansion, simplifying 
insofer as possible. This means that the rules are essentially multipli- 
cative in nature. The process of simplification tends to conceal this 
fact, so that they appear to be at least partially additive. 


The multiplicative nature of the rules is apparent from the view- 
point of conformal mapping. Thin-airfoil theory gives a formal series | 
for the mapping of the profile onto a straight line, which fails near a 
round nose. With Lighthill's rule, thin-airfoil theory is relied upon 
only to map the profile onto an initially osculating parabola. The radi- 
cal in equation (3b) is the speed on the parabola; the remaining factors 
represent the speed ratio associated with this intermediate mapping, by 
which it is to be multiplied. 


From another point of view, Lighthill's rule implies that the airfoil 
is represented by sources distributed over the surface of a parabola, 
rather than along a straight line. To second order, the perturbations 
induced by the sources are proportional to the local speed of the basic 
flow past a parabola, which indicates again the multiplicative nature of 
the rules. This viewpoint remains valid in the subsequent extensions to 
subsonic and axially symmetrie flows, where conformal mapping is inap- 
plicable. 


ROUND-NOSED ATRFOILS IN SUBSONIC FLOW 


The preceding results can be extended to compressible flow by using 
solutions for subsonic flow past parabolas and ellipses. Although no 


exact solutions are known, it seems that adequate approximations can be 
found so long as the flow is purely subsonic. 
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At stagnation points in subsonic flow, thin-airfoil theory predicts 
infinities of the same order as in incompressible flow. Although these 
infinities can be avoided by using the exact condition of tangent flow 
at the surface (which in incompressible flow gives the true solution), the 
result is nevertheless incorrect near stagnation points. Thus, Hantzsche 
and Wendt (ref. 10) work with the perturbation stream function and impose 
the exact tangency condition by means of conformal mapping, but even their 
first approximation gives not zero speed at stagnation points but a veloc- 
ity u/U = -M?(1 - M“)"1/2, This is primarily a defect associated with 
the stream function, but even using the velocity potential in the same 
way gives results which are no more than qualitatively correct near stag- 
nation edges. 


An alternative procedure for treating subsonic flows by successive 
approximations is the Janzen-Rayleigh expansion in powers of M°. The 
two methods are known to complement each other, being, in fact, simply 
different series representations of the true solution. In particular, 
the Janzen-Rayleigh result is most accurate near stagnation points, where 
thin-airfoil theory breaks down. It is therefore suited for the present 


purpose. 


Consider again an uncambered round-nosed airfoil at zero angle of 
attack. Probably the second-order thin~-airfoil solution is the most that 
will be treated in practice; then it is enough to consider subsonic flow 
along the axis of a parabola (sketch (a)). Imai has recently calculated 
(ref. 11) the Janzen-Rayleigh solution for a parabola including terms 
in Mf, and has given the surface speed in essentially the form 


F = fo + Mf. + (y + 1)Mffa + Mfg + O(M) (13) 


Here the f, are increasingly complicated functions of x/p (where p 
is the nose radius), fo being the incompressible result of equation (1). 


Far behind the nose, Imai's result reduces to 





: M4 rmt) + O(MÎ, p?/x2) (14a) 


This should agree to this order with thin-airfoil theory and indeed it 
does; second-order thin-airfoil theory (ref. 12) gives 





"a2 P o(y + 1)M2 + bB?) 2, 2 
= = 1 is L+M na |+ otee (14) 
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and. the reduced result of Imai, equation (lka), is just the curtailed 
expansion of this in powers of M°. Comparing these two expressions | 
suggests a modification of the Janzen-Rayleigh solution that will include 
the full second-order thin-airfoil result: | 


= ENE x) = fo + ME f +(+ De 34 x fe + M*(fg - fx) + O(M®) (15a) 


where, according to Imai's solution,* ` 


1/2 
Ta = let = Sin 2 where 0 = cos”? > 
j x + 9/2 2? o P + 2x 


fy = -= 2 (cos 2) E tan ons 2 sin 9 log (2 cos 2) + @ cos J 
l 2 2 p 2 


: 5 | | 
E cc) ca fy dee _ 
fs = 33 ( cos 2) |x tan 5 13 sin 6. Wsin @ - 3 sin 29) 


log (2 cos E) + 20(cos 8 - 3 cos 2) | 


e-n- -p(o 8) |-GE- 5) amore 
a “i1 = cos es eg ae sin 9 + Z- sin 20 + 7 sin 20 
(15b) 
108 (2 cos 2) - Zis - 12 cos 6 + 7 cos 20) + 


2 f 
(3 sin @ + 8 sin 20 +- 3 sin 30) + 

L 2 9 
5isin 68 =- 4 sin 26 - 3 sin 36) log (2 cos e fi 


H5 cos @ + 8 cos 20 + 3 cos 30)10g( 2 cos 3 z 


o 
2(2 cos @ + cos a 108 ( 2 cos 3) 29 | 


expression for fis the correction has been confirmed by Prof. Imai. 
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This mođification seems to be a natural one; for example, the combination 
(fa - fı) is simpler than its original counterpart fa. A short table 

of the function Q(x/p,M) calculated from this approximation? is given 
below for y= 7/5. 





Applying the argument used before to equations (14) and (15) shows 
that the second-order thin-airfoil solution for any uncambered analytic 
profile of nose radius 9 in subsonic flow at zero angle of attack is 
converted into a uniformly valid second approximation by the rule 


d2 _ Mas „O + 1)M + hp? 
U (Zsu) E- aaa 2 | (16a) 


Here Q(x/9,M) is the exact subsonic solution for a parabola. It can, in 
principle be found to any desired accuracy since the Janzen-Rayleigh 
method is believed to converge for purely subsonic flows and, for axial 
flow past a parabola, the flow is purely subsonic for any free-stream 
Mach number less than unity. For practical purposes, the M* approxima- 
tion of equation (15) that is tabulated above will probably suffice. For 
application to linearized theory, this rule reduces simply to 


SGE on 


Which yields a uniformly valid first approximation. 





This last rule actually applies as a first approximation to any 
round-nosed airfoil at the ideal angle of attack. Corresponding rules 


“These values have been verified by comparison with a more extensive 
table of the coefficients OCENETSRE in equation (13) that was kindly 
communicated by Prof. Imai. 

SA convergence proof was reported in reference 13, but has never 
been published. | 
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for lifting and cambered airfoils at other angles of attack can be derived 
from the Janzen-Rayleigh solution for a parabola with circulation.’ This 
can be extracted to order M? from existing solutions for inclined 
ellipses; however, these solutions have not yet been carried as far as 

the M* terms, which are the first to involve y and are probably 
important for numerical accuracy. 


The pressure coefficient must be calculated from the velocities 
according to the full Bernoulli equation for isentropic flow: 


; o 


Y 
— 
y -1 


2 a j | gS 
SHARP -NOSED AIRFOILS 


Aside from round noses, sharp edges with finite vertex angle arise 
most commonly in practice. At sharp edges, thin-airfoil theory predicts 
logarithmically infinite velocities in place of stagnation points and, 
for the circulatory component, predicts a leading-edge singularity that 
exceeds its actual strength. From a practical point of view, it is 
probably worthwhile to correct these shortcomings only at sharp leading 
edges because at trailing edges the details will ordinarily be masked 
by viscous effects. In any case, the correction is significant over a 
much smaller region, for a sharp edge than for a round one. 


Consider first incompressible flow past a sharp-nosed symmetrical 
airfoil at zero angle of attack (sketch (e)). Clearly the flow in the 


yo ae immediate vicinity of the nose is 
/ like the flow in an angle.8 By con- 
/ Ko formal mapping, or separation of 
LOTT ži variables in polar coordinates, the 
p \ , surface speed is found to be given 
\ y 
N / 5 
———— q pa 
Sketch (e) = = ox” (18) 


TBecause of the nonlinear nature of compressible flow, here the 
"ideal" and circulatory components can no longer be treated separately, 
Which somewhat complicates the rules. Even in incompressible flow, if the 
two components are treated together, there results instead of equation (3b) 
an alternative rule that is more complicated but, of course, entirely equi- 
valent up to terms of second order. In particular, this alternative rule 
requires that the coefficient k of equations (4) and (5) be known. Hence, 
in the rule for subsonic flow at other angles of attack than the ideal, it 
will be necessary to determine k as the limit of the coefficient of x71/2 
in the thin-airfoil solution as the airfoil thickness ratio vanishes 
(corresponding to transition to a true parabola near thé nose). 


&See reference 14 for a rigorous discussion of this point. 
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where 5 is the semivertex angle. In fixing the constant c, the diffi- 
culty that in the angle flow the velocity increases indefinitely upstream 
can be circumvented by requiring that at any point the velocity must 
approach that of the free stream as the angle 6 tends to zero. Thus, 
it is seen that ce is unity except for terms of order ô. 


The connection with thin-airfoil theory follows from the fact that 
for small WV 
xV = eY 108 X 214 Vlog x+ Sv? log@ x+... (19) 


though not uniformly near x = 0. Hence, the thin-airfoil series for 
flow in an angle is, to second order, 


2 
= eft + È tog x + (tog x + È 10g? x) (20) 
K 2 


Gha 
j 


and it is clear how the spurious logarithmic singularities arise. 


Comparing these expressions gives a rule that renders the second- 
order thin-airfoil solution for any sharp-nosed profile uniformly valid: 


= 5 Tn 
q> Pa = de 5 "qt 52 l 2 
T = x“ | y ~ x U log. x= -2 log x = 5 log x (21a) 


The corresponding first-order rule is 


aa.: a= dog x) (21b) 


These rules, unlike those for round noses, can be extended to indefinitely 
high order simply by retaining more terms in equation (20). The reason 
for this is that the region of nonuniformity is exponentially small, so 
that the shape of the profile (if it is composed of analytic arcs) enters 
only through its initial angle. 








The second of these rules applies 
as a first approximation to any sharp- 
nosed airfoil at its ideal angle of 
attack. At other angles the flow 
includes a circulatory component, which 
behaves near the edge like that around 
the angle shown in sketch (f). For this 
component, conformal mapping gives the 
surface speed as 





Sketch (f) 
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T=-285 
=k sin dg x 2(t-5) e +h 8ing +2 log x +. bi. i >) (22) 


tx 


The constant k, which need not be identified, is of order unity, so 
that logarithmic terms arise in thin-airfoil theory only in the terms of 
second order in a, and 5. The rule for correcting this circulatory 
component is, to first order, 





GI] a 


5 
GM _ ,e(n-5) To 1 
7 = = (23a) 
and to second order 
a LA ee $ i 1? 
de (x-8 də = age. 


Note that according to thin-airfoil theory, the circulatory component 
involves speeds with an inverse square-root singularity, the actual sin- 
gularity being slightly weaker. Here the rules, in contrast to those for 
round noses, are different for the two components. Hence, the thin- 
airfoil solution must be treated by splitting off the terms which are 
singular at least as x12, applying equation (23) to this circulatory 
component, and then applying equation (21) to the remainder. 


It is clear that corners on a profile elsewhere than at the leading 
or trailing edges will lead to logarithmic singularities in thin-airfoil 
theory that can, for incompressible LOW; be treated by a modification of 
this procedure. 


At subsonic speeds, the flow at the ideal angle of attack could be 
treated by considering the Janzen-Rayleigh solution for flow in an angle 
(which has not yet been calculated). However, at other angles the Janzen- 
ities that are tolerable only in an incompressible fluid. The present 
comparison method could be applied only if a erly transonic solution 
were known for an inclined wedge. a 


COMBINED EDGES; AN EXAMPLE 


The rules derived above for a single edge at the origin can clearly 
be combined to treat leading. and trailing edges together. In order to 
use the previous results directly it is essential, as Lighthill has 
pointed out, that the x axis be chosen to pass through both edges. It 
is also essential to bear in mind the multiplicative nature of the rules. 
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Thus, consider incompressible flow past an airfoil with round lead- 
ing and trailing edges of radii 9, and Pp located (without loss of gen- 
erality) at x = -l and x = 1, respectively. Applying equation (3a) 
twice, identifying x successively with (x + 1) and (1 - x), and then 
simplifying to keep no more than second-order terms (recalling that the 
p are of second order) gives 








E 1/2 
a Lex (< 7. wo I f 
U U hie, bie Gi 


L 
l =- x? +4 5 Pall- x) + =pp(l+ x) 


(Note that the simplification has included dropping a term PaPp/} in 

the denominator of the radical, since its contribution can be shown to be 
of higher order.) This rule converts the second-order thin-airfoil solu- 
tion into a uniformly valid second approximation (in the sense discussed 
previously) for airfoils whose initial and final camber is zero, but gives 
only ea first approximation near a cambered edge. 


In the same way the previous rules can be combined for higher approx- 
imations and for any combination of round and sharp edges. 


As an application of the rule for combined round edges, consider 
incompressible flow past an ellipse of thickness ratio € at zero angle 
of attack. The formal second-order solution for surface speed is 


"a2 1 x2 
ge re eS ee (25) 


= 1/2 
Se (ais ( E 2) 
7 te) l+re+se (26) 


and it is satisfying to see that aside from the radical, the expression 
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has been simplified. This result is compared in sketch (g) with the 


L5 





—— Exact 


eseese Linearized theory 


---= Second-order theory 
© Corrected second-order theory 


90 
8, deg 


Sketch (g) 


exact solution and with the formal first- and second-order solutions for 
an ellipse of thickness ratio e= 3/10. 


APPLICATION TO FINITE WINGS 


The edges of a finite wing are described as subsonic or supersonic 
according as, the component of free-stream velocity normal to them is 
subsonic or supersonic. At subsonic edges thin-wing theory breaks down 
just as it does for..two-dimensional airfoils in subsonic flow.9 These 
defects can be corrected by applying the preceding rules to the components 
of velocity in the plane normal to the edge; for, as noted by R. T. Jones 
(ref. 2), for sufficiently thin wings, the flow field in the immediate 
vicinity of the edge is cylindrical (if the plan form and edge radius are 
continuous) according to the simple Sweepback concept. Since thin-wing 


° Thin-wing theory fails also at round supersonic edges. This defect 
could probably be corrected by comparing with the solution for supersonic 
flow past a parabola, which has not yet been calculated. 


3D 


w’ 
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theory is a series expansion in powers of thickness, its coefficients 
can be evaluated at vanishing thickness. Hence, the sweepback principle 
applies here to any order. i 


As an example, consider supersonic flow past a flat elliptic cone 
lying inside the free-stream Mach cone. The Linearized thin-wing solution, 
first given by Squire (ref. 15), predicts constant pressure over the entire 
surface. Hurley has attempted to improve this result near the leading 
edges by adding, along the foci of the ellipse, conical line sources whose 
strength is determined by imposing the exact tangency condition at the 
leading edge (ref. 16). He has verified that the resulting shape is 
approximately the desired one, though somewhat sharper near the leading 
edge. For the thickest wing that.he considered, Hurley's pressure dis~ 
tribution is compared in sketch (h) with the first-order result of the 
present method. The present solution is presumably more accurate since 
the shape is not distorted and compressibility effects are accounted for 
more exactly. 


4 


--+- Linearized theory 
` —— Hurley (ref. 16) 
tan’ 1S —— Present method 


Cp .2 
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Sketch (h) 


BODIES OF REVOLUTION 


For bodies of revolution, the counterpart of thin-wing theory is 
the slender-body theory introduced by Munk for the crossflow due to angle 
of attack, and later extended to cover also the flow at zero angle. For 
subsonic speeds this theory, too, fails near stagnation points, where it 
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predicts infinite rather than zero speeds. Furthermore, for the crossflow 
it predicts a nonexistent singularity at round ends and fails to predict 
an actual singularity at sharp ends. It can be corrected by simple rules 
quite analogous to those already given for wings. Only incompressible 
flow will be considered here, the extension to subsonic speeds requiring 
the solution for simple nose shapes according to the Janzen-Rayleigh 
method or otherwise. : 


Any round-nosed body of revolution can be approximated near its nose 
by a paraboloid of revolution. Munk's rule gives the remarkable result 
that the surface speeds for axial flow are identical on parabolas and 
paraboloids of revolution, and that the additional velocities due to angle 
of attack are proportional. Hence, Lighthill's rule for round=-nosed air- 
foils, equation (3b), applies also to bodies of revolution. It yields, 
however, only a first approximation because the disturbances, being 
smaller than for airfoils, are of the order of the nose radius (i.e., of 
order ¢2) rather than of order e€. A second-order rule can be obtained 
by fitting the nose more accurately with an ellipsoid of revolution. 
Again, Munk's rule shows that the speed distributions are proportional 
on ellipses and ellipsoids of revolution of the same thickness ratio. 
Hence, the fourth-order rule for airfoils, equation (11), is the second- 
order rule for round-nosed bodies of revolution. 


For sharp-nosed bodies, the nature of the flow near the tip can be 
found by introducing spherical polar coordinates (sketch (1)). Separating 
variables leads to a solution in 
terms of Legendre functions for the 
axial flow and associated Legendre 
functions for the crossflow, the 
functions being of nonintegral order. 
Requiring, as in plane flow, that 
the solution reduce to a uniform 
stream as the cone angle © vanishes 
shows that on the surface in the 
immediate vicinity of the tip the 





Sketch (1) radial velocity component has the form 
= = c cos a xf" + k sin a cos w xB7t (27a) 


and the azimuthal component, the form l l n 


q , ` 
y = -k sina sin w x0" 


(27b) 


Here, again, c is a constant that is unity to first order and k isa 
constant: of order unity. The tangency condition requires that m and n 
be the first roots of 


bey 


4 . 


|. pal 
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|! 
O 


P,,’ [cos(x - 5)] 
(28a) 
Ppt’ [cos(x - 8)] 


It 
© 


where the derivatives of the Legendre functions are considered as func- 
tions of mand n. Methods for calculating these roots numerically have 
been discussed by Pal (ref. 17) and by Siegel et al. (ref. 18). Using 
a formula due to Schelkunoff (ref. 19), it can be shown that to a first 
approximation for small cone angle 5 


m= 1 +5 8 + 0(8* log d) 
(28b) 
1. 
n= 


= 1-58? + 0(8° log 8) 


Since n ~- 1 is slightly negative, the velocities associated with the 
crossflow are only weakly infinite at the tip, in contrast to the nearly 
inverse square-root singularity associated with angle of attack for a 
sharp-nosed airfoil. 


Expanding equation (27a) according to equation (19) leads again to 
rules for correcting slender-body theory. As in plane flow, the axial 
Flow and crossflow must be treated separately. The rule for correcting 
the axial component is found to be, to first order, 





qı = qi L 
2. (3 - F 88 log x ) (29a) 


This could be extended to second order by determining the next term in the 
series for m. The rule for correcting either the radial or azimuthal 
component of velocity associated with the crossflow is, to second order, 


qo n=1 ( "ds l eo l ) 
Se 2 sis 
7 =X es log x (29b) 


The corresponding first-order rule is the same without the logarithmic 
term. 
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ALTERNATIVE FORMS OF RULES FOR ROUND NOSES 


1/2 
The radical (is appearing in the rules for round noses in 
incompressible flow is, for a parabola or paraboloid, the cosine of the 
surface angle ņ measured from the x axis. Therefore, the first-order 
rule for airfoils, equation (8), can be written alternatively as 
= 1 S $ 
B = cos n (E) > 
1 + tan? n 
and this is Riegels! rule (ref. 9). Aside from its simplicity, it has an 
unexpected advantage over Lighthill's form. _The two rules differ only by 
higher-order terms for shapes other than a parabola, so that there is no 
a priori reason to prefer either. However, it happens that Riegels! form 
gives the exact result for an ellipse as well as for a parabola,?° whereas 
Lighthill's does not. Since airfoil profiles ordinarily resemble ellipses 
more than parabolas, Riegels' form might be expected to be more accurate, 
and this appears to be the case. For example, sketch (j) shows that 
Riegels! rule gives a much more accurate first-order solution than 
Lighthill's for an NACA 0012 airfoil at zero angle of attack. 
LZ 


(30) 





gq 


L2 NACA 0012 


Exact (ref. 8) 


Linearized theory l 
First order, Lighthill's rule 
First order, Riegels" rule 


40 


0 2 4 6 a LO 
Fraction of chord 
Sketch (j) 


This is clear from Munk's rule (together with the fact that line- 
arized theory predicts a constant speed on an ellipse that is equal to 
the true maximum value). 
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Despite its simplicity and remarkable accuracy, however, Riegels! 
rule must be regarded as less basic than Lighthill's from a theoretical: 
point of view. Only Lighthill's rule can, in principle, be extended to 
indefinitely high order, to subsonic speeds, or to other nose shapes. 

A further defect of Riegels' rule is that the nose exerts a spurious 
influence even at remote points if the local airfoil slope is appreciable. 
To this extent it fails to render the solution uniformly valid. 


Nevertheless, the practical advantages of Riegels! rule are so great 
as to suggest artificially modifying the higher-order rules accordingly. 
Thus, Lighthill's second-order rulé for uncambered airfoils, equation (3b), 
can be rewritten by analogy with Riegels' rule as 





= 17 tf 
2 de L 
-J 7 cos a( v + = tan? n) (31) 


Both forms are so accurate that graphical comparison is scarcely possible. 
Instead, they are compared in the following table with the exact values 
for a 14-percent-thick symmetrical Joukowski airfoil at zero angle of 
attack: 


Exact 
(ref. 20) 





The superiority of the alternative form is now confined to the immediate 
vicinity of the nose (which has a radius of 0.023 chord). Similarly, for 
cambered airfoils, the second-order rule can be recast in the form 





om Ta 
do 2 
a = cos() - k) | 7 


Where Kk is the initial angle of the camber line; and comparison with 
its counterpart in equation (10) illustrates the greater simplicity of 
the alternative forms, which is perhaps their chief virtue in the higher- 
order approximations.’ 


+> tan®(q - e) (32) 


Again, both these rules are exact for ellipses, whereas their pre- 
decessors are not. Since this is so, the fourth-order rule (which results 
from considering ellipses) is obtained simply by retaining another term 
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in the formal series expansion of the cosine. For initially uncambered 
airfoils this leads to 








nn "q 1 "a. 1 
Ga _ 4 ,i fə EEEE: eee 
g = cos 7 ( yt ST tan" y -g tan n) (33) 


and it can be verified that up to terms of fourth order this agrees with 
the less compact form of equation (11). The corresponding rule for 
ecambered airfoils is 


re Wi tt 2 
4 
— = cos() - 








= = = tan=(n - k) - = tan*(n - e) (34) 


For round-nosed bodies of revolution, equation (31) is the first- 
order rule and equation (33) the second-order rule. These are exact for 
ellipsoids of revolution as well as paraboloids. It can therefore be 
anticipated that they will be generally more accurate than the previous 
forms. 


The rule for subsonic speeds can also be artificially manipulated 
in an attempt to preserve the advantages of Riegels! rule. Thus, equa- 
tion (16a) is replaced by 








do ~ "g" 1 M2 LM 
2 | 7 pOf H) tan? n (35) 


where Q(n) is the subsonic speed.ratio on a parabola expressed in terms 
of the surface angle yn. It is related to the function Q(x/p) of egua- 
tions (15) b 


An) = aG cot? a) (36) 


because for a parabola p/x = 2 tan? ne This rule differs from its pre- 
decessor only by higher-order terms and is designed to reduce to equa- 
tion (31) at zero Mach number. 


Ames Aeronautical Laboratory 
National Advisory Committee for Aeronautics 
Moffett Field, Calif., Aug. 16, 1954 
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APPENDIX 
SYMBOLS 
pressure coefficient (pressure increment divided by free- 
stream dynamic pressure) 
initial slope of camber line divided by e€ 
coefficient that is unity to first order 
coefficients in expansion of airfoil thickness 
coefficients of Janzen-Rayleigh expansion for Q 
coefficient of order unity 
free-stream Mach number 


first root of Py'[cos(x - 8)] = 0, where Py,(cos 0) is 
the Legendre function 


first root of Pni*[eos(xn - 5)] = 0, where Pmè(cos @) is 
the associated Legendre function 


speed a on parabola in subsonic axial flow as function 
of x/P 


speed ratio on parabola in subsonic axial flow as function 
of surface angle n 


speed at body surface 

uniformly valid approximation to q 

formal thin-wing or slender~body solution for q 
first-order, second-order, etc., approximation to q 
free-stream speed 

component of surface velocity in x direction 


coordinate along airfoil chord or body axis, usually 
measured from edge 


coordinate of airfoil surface perpendicular to chord 
angle of attack 


angle of attack measured from ideal angle 
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B (1 - M?)7/* 

yY adiabatic exponent of gas 

© semivertex angle of pointed airfoil or body 
E thickness parameter of airfoil 

7 surface angle measured from x axis 

9 polar angle 

K initial angle of camber line 

A (see eq. (11)) 

p edge radius of round-ended airfoil or body 


W azimuthal angle 
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